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Abstract. We consider the influence of a power-law deviation from the critical 
coupling such that the system is critical at its surface. We develop a scaling theory 
showing that such a perturbation introduces a new length scale which governs the 
scaling behaviour of the density profiles as well as the finite-size behaviour of the 
surface properties. Exact results are obtained for the Ising quantum chain when 
the perturbation varies linearly whereas the quadratic perturbation is mainly studied 
numerically. The scaling theory is well confirmed in both cases. 



1. Introduction 

Real systems are in general inhomogeneous due to the presence of impurities, line defects, 
boundaries, aperiodicity or disorder. Such inhomogeneties may have a strong influence 
on the critical properties in the vicinity of a second order phase transition [H El El HI [5] . 
Depending on the relevance of the perturbation introduced by the inhomogeneity, the 
universality class governing the critical behaviour may change; in some cases, the critical 
singularities may even be suppressed. 

Besides this, the presence of spatially varying external fields (magnetic, 
gravitational or thermal) also influences the critical properties of a system. The effect 
of gravity on phase coexistence was studied in [61 IE] • Phase separation induced by a 
thermal gradient was considered in [9], [TDl [TTj . The application of a gradient perturbation 
was used as a tool for high-precision estimates of percolation threshold and exponents 

in nana nana. 

In general the effect of a spatially varying field on a critical system is to smooth 
out the critical singularities. The reason for this is that the perturbation generated by 
the inhomogeneous field leads to a departure from the critical point which introduces a 
finite length scale in the problem [32] . One may think about the usual finite-size scaling 
theory as a particular example where the spatially varying field is a confining potential 
of a box-like type [T6j . 
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For a field generating a power-law deviation from the critical point, A(z) = gz w , 
the typical length scale introduced around the critical locus, z = 0, is £ ~ \g^- u /( 1 + l ' u> ) 
where v is the correlation exponent. Far away from this critical region, for \z\ ^> £, 
the correlation length is small compared to the scale of the spatial change and the 
local properties at z are determined by the corresponding homogeneous system with 
a constant deviation to the critical point given by the value A(z). Based on these 
observations, a scaling theory for the density profiles in the presence of gradient-field 
inhomogeneities was developed in [T7] . The validity of the scaling theory was checked at 
the mean-field level within Ginsburg-Landau theory and by comparing its predictions 
to the exact solution of an Ising quantum chain in a linearly varying transverse field. 

The same scaling predictions were tested for the Ising quantum chain with a 
linear-to-flat profile of the transverse field [18], for a spin-1 ferromagnetic Bose-Einstein 
condensate in a space-varying magnetic field [19] and for a classical Ising lattice gas 
confined in a two-dimensional harmonic trap |20j. The entanglement entropy for 
Ising and XX quantum chains with linearly varying fields or couplings was studied in 
[2 lj . Quite recently the trap-size scaling was examined for the XY and Bose- Hubbard 
quantum chains [22]. 

One may mention a formal analogy noticed in [18] between the spatial 
inhomogeneous case discussed here and the quench dynamics situation where a system 
is driven through a critical point at finite rate. During the quench, far away from the 
critical point, the system is able to adjust its state adiabatically as far as its relaxation 
time is much shorter than the typical time associated with the quench. On approaching 
the critical point, due to the critical slowing down, the typical quench time becomes 
much shorter than the system relaxation time and the evolution switch approximately 
to an impulse regime. Close to the critical point, the typical time scale during which 
the state is frozen (impulse evolution) is given by equating the system relaxation time 
(inverse of the first gap for a closed quantum system) to the typical quench time scale. In 
the quantum case this leads to the self consistent equation r = | l/e(r) \ uz where e(t) is the 
deviation to the critical point at time t, z the dynamical exponent and v the correlation 
length exponent. For a power law deviation e(i) = gt u one obtains r ~ g-zv/(i+wxv) 
in complete analogy with the spatial case. Drawing on this analogy, we see that to 
the adiabatic evolution regime corresponds a "super-static" regime where, locally, the 
physical properties are determined solely invoking a local equilibrium hypothesis. On 
the contrary, on approaching (either in time or in space) the critical locus, the local 
equilibrium hypothesis (time- or space-like) breaks down due to the divergence of the 
typical time (relaxation) or space scale (correlation length). 

All these studies were made in reference to the bulk critical behaviour, i.e., when 
the locus of the critical field value is deep into the bulk of the system. However, if 
one shifts the critical locus toward a free surface, when close enough to it, the scaling 
behaviour of local quantities should depend on the surface critical properties. In the 
present work we consider such a situation for a semi-infinite system with an ordinary 
surface transition. The paper is organized in the following way: in the next section 
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we present the scaling arguments for a semi-infinite system with a power-law deviation 
from the critical point when the critical locus is at the surface. The third section is 
devoted to the semi-infinite Ising chain in a linearly or quadratically varying transverse 
field. In the linear case we solve analytically the problem while in the quadratic case 
the results are obtained through an exact numerical diagonalization. In both cases our 
findings are in perfect agreement with the scaling predictions. Finally we summarize 
our results and indicate some possible developments in the last section. 



2. Scaling arguments 

In the next sections we study either linear (gradient) perturbations, for which exact 
results can be obtained, or quadratic ones. Thus here we shall examine the scaling 
behaviour for an arbitrary power law variation of the couplings. 

The system is semi-infinite, occupies the half-space z > with free boundary 
conditions at z = 0. The perturbed coupling deviates from its critical value K c a^§| 

K{z) - K c = A(z) = gz" , uj > (2.1) 

Thus the coupling is critical at the surface z — 0. The system is ordered (disordered) 
in the bulk when g > (g < 0). The perturbation introduces a crossover length 
I giving the width of the surface region which diverges when g vanishes. With a 
decaying perturbation (uj < 0) the discussion of the scaling behaviour is similar but 
the physics is quite different since the bulk remains critical while the surface critical 
behaviour is mofified when the perturbation is relevant (0 > uj > — 1/V) or marginal 

(u = -i/i/)PHl E1I2SI ESI E3E2B]. 

Changing the length scale by a factor b, the thermal perturbation A(z) with scaling 
dimension y t — \jv where v is the correlation length exponent, transforms ajjj] 



g'^ = b^gz" = g> ft j , (2.2) 



so that 



9 



' = b w+1/u g . (2.3) 



As mentioned above g grows under rescaling and the perturbation is relevant when 
uj > —1/v. 

The crossover length transforms as 

t = i = i{g') = t (V+^g) , (2.4) 

with b = \g\~ u /( 1+UUJ \ one finally obtains 

£ = A\g\-"^ 1+vu \ (2.5) 

§ For g < equation (|2.ip leads to negative values of the coupling for sufficiently large z. The value 
of g and the size of the system should be chosen in order to avoid the occurence of negative values. 
|| The following discussion extends to other types of perturbation provided the appropriate scaling 
dimension replaces yt- 
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for the crossover length introduced by the thermal perturbation. The result of equation 
(2.5) in [T7] is recovered when u — 1. 

As in [IT] this result can be obtained self-consistently p2J El] by relating £ to the 
local value of the correlation length £(£) so that 

e<xZW<x\A(l)\-"<x(\g\n- v , (2-6) 

which gives back ( 12.51) when solved for £. 

Let us nown examine the scaling behaviour of the profile (p(z,g) associated with 
some bulk operator ip with scaling dimension x v . This may be the order parameter m 
or the singular part of the energy density e. It transforms as 

<pf = <p(z , ,g r )=b B *<p(z,g), (2.7) 
so that, using (12.31) 

<p(z,g) = b-*><p(l i P +1 '''g) . (2.8) 

With b = \ g \-"/( 1 + vu ) oc £ one obtains 

<p(z,g) = (\g\ v l^ v ^ z ) oc £- x ^(A^) , (2.9) 



where the last expression follows from (12. 5p . The prefactor exhibits the finite-size 
behaviour expected for a critical system with size £. In the same way, for a surface 
operator ip s with scaling dimension x Vs at the ordinary surface transition, in a system 
with size L, one obtains: 

(f s {L,g) = |^W(i+M$ s (|^/(i+^) jL ) ocr 2 ^$ s ^4^J . (2.10) 

When L £, the effect of the perturbation is negligible and the surface operator displays 
the usual finite-size behaviour for a critical system ip s (L,g) ~ L~ Xvs . As a consequence, 
the scaling function behaves as: 

$ s (w)~w-^, u<l. (2.11) 

When L/£ — > oo, <p s no longer depends on L, so that: 

lim $ s (w) = const , lim ^ s (L,g) oc r x * = const \g\ vx ^/a+^) . ( 2 .12) 

Such a behaviour is expected with e s for any sign of g and with m s for g > 0. 

In the surface region z £ the profile in (12.91) . when non- vanishing, should scale 

as £~ x ^ 

<p(z, g) oc e- x *$(A |) oc r a « , | < 1 , (2.13) 

leading to: 

~ ; w < 1 . (2.14) 

The order parameter m(z,g), with scaling dimension x m = /3/u, is non- vanishing 
in the ordered phase g > 0. Then, outside the surface region when z>f, one expects 
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the same local critical behaviour as for the homogeneous system with the same value of 
the coupling, i.e., 

m(z, g) oc A(zf oc (gz^f^ , g > . (2.15) 

for not too large values of A(z). Comparing (12.151) to the first expression of the profile 
in (12.91) with (p = m, one obtains the asymptotic behaviour of the scaling function: 

$ m (u) ~ u^*™ , u > 1 , g > . (2.16) 



3. Quantum Ising chain in a varying transverse field 

3.1. Quantum Hamiltonian 

We consider the quantum Ising chain in a varying transverse field with Hamiltonian [29] 

L-l L 

i=i i=i 

where the a x,z are Pauli spin operators. This Hamiltonian governs the extreme 
anisotropic limit of the classical two-dimensional Ising model on the square lattice 
[30| |3"H 32J. As explained in [17J, when the classical system has horizontal couplings 
K2 and varying dual vertical couplings K*(l) = K^(l — gl w ), the row-to-row transfer 
operator T takes the form 

T = 1-2K 2 H. (3.2) 

when K\ — > 0, K2 — > while /i = K\jK 2 remains constant. Universality holds in 
this limit and the critical properties of the quantum chain are the same as those of the 
two-dimensional classical system with pQ: 

v = 1 , x e — 1 , x m = 1/8 , x Gs = 2 , x ms = 1/2 (ordinary transition) . (3.3) 

In the following we take h — 1 and the quantum chain is truly critical when g = 
and L — >• 00. 



Diagonalization 

Under a Jordan- Wigner transformation [33] the Hamiltonian (13.11) is changed into a 
quadratic form in fermion creation and destruction operators which is diagonalized 
through a Bogoliubov transformation [3H [32] . Then 

n = Y,e q U Vq -l-) , (3.4) 

where rj^ (r) q ) are diagonal fermion creation (anihilation) operators. The fermionic 
excitation energies e q are obtained through the solution of one of the following eigenvalue 
equations 

<f>g{l ~ 1) + [e 2 q - 1 - hf\ (j) q {l) + hi <f>g{l + 1) = , 

hi %{l - 1) + [e\ - 1 - /if] ^,(Z) + h l+1 if, q (l + 1) = (3.5) 
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0. 

= 



with the boundary conditions 

(? (O)-0,(1) = O, <f> q (L + l) 
V 9 (0) = 0, %{L)-iP q (L+l) 

In these equations ip q (l) and (j) q {l) are the normalized eigenvectors 
The scaling variable of the quantum chain with v — 1 is 

|g|!V(HW>; = |^|l/(w+l)^ _ 



(3.6) 



u 



In the scaling limit g — > 0, L — > oo while keeping # 



«(/+!) -u(l) 



\9\ 



and w(L) = \ g \-UM^)]g^ 



(3.7) 

constant, we have 
► oo. Thus the system is 



ii/(«+i) 

semi-infinite in u which is a continuous variable. The finite-difference equations in (13.51) 
can be rewritten as differential equations in the reduced variable u. Keeping the leading 
contributions in g, of order \g\ 2 '^ +1 ', one obtains the following differential equations: 



du 2 

d 2 V> 
du 2 



+ 



+ 



ii/V+i) 



e 



\9 



sign(#) uju^ 1 
- sign(g) uju^' 1 



u 







d0 
du 



0. 



oo 



u=0 



u 



2^ 



ip(u)=0, </>(0) 



0, ** 
du 







0. (3.8) 



When g changes sign the two equations are exchanged but the boundary conditions 
remain the same. 



4> s>Vs . 

4>i»Yi • 
i-¥ • 



(j) 5)¥5 . 

<t>3>Y 3 • 
4> 2 .V 2 ■ 

♦i.¥i ■ 



g<0 



<t>0 



g>0 



Figure 1. Excitation spectra of the quantum spin chain. When g > (ordered phase) 
the lowest excitation vanishes and the ground state is degenerate. 



When oo = 1 we have: 



d 2 

du 2 


-Isl 


d 2 ^ ^ 
dw 2 


~e 2 


.\g\ 



- sign(p) - m 2 





' du 



0, 



00 



«=o 



V»(«) = o, v(o) = o. 



d^ 
dw 



0, (3.9) 



with u = g 1 / 2 ^ where z > is the continuous coordinate along the quantum chain in 
the scaling limit. 
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These equations are eigenvalue equations of the harmonic oscillator forrrjjj] 



du 2 



+ (2n + 1 -u 2 )xn(u) = 



(3.10) 



with eigenf unctions x n (g 1//2 z), vanishing at infinity and normalized on z > 0, which are 
given by 



Xn{U) 



/2 H n (u) 



HJu) 



'-l) n e u 



(3.11) 



2 n n\ \nJ uy " ' lx ' v ' du n 

where H n (u) is the Hermite polynomial of order n. In ( 13.91) . at u — 0, the Neumann 
boundary condition for <p selects even values of n whereas the Dirichlet boundary 
conditon for ip selects odd values of n. Thus for the ordered system (g > 0) one obtains: 

£ P =V^P9, M u ) = X2 P (u), ip p (u) = X2 P -i(u) , p= 1,2,3, ■■■ 

e = 0, (j) (u) = Xo(u) , tp (u) = 0. (3.12) 



The vanishing excitation eo is a consequence of the twofold degeneracy of the ordered 
ground-state. The behaviour is different for the disordered system (g < 0) with a non- 
degenerate ground state. Then the lowest excitation is non-vanishing: 

e p =y/2(2p+l)\g\ 



<f>p(u) = X2p(u) , 1p p (u) = X2p+l(u) , p 

The lowest part of the excitation spectra is shown in figure [TJ 
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Figure 2. Scaling of the magnetization profile for ui — 1 (left) and uo — 2 (right), the 
lines indicate the slopes expected from (|2.14p for z <g; I and (|2.16p for z>f. 



3.3. Magnetization profile and surface magnetization 

When g > the Z2 symmetry is spontaneously broken in the thermodynamic limit. In 
a finite system it can be broken by fixing the last spin (<r£ = ±1 when hi = 0). Then 
the magnetization is non- vanishing. At 11, it is given by the off-diagonal matrix element 
|((j|af|0)| where |0) is the ground-state of H and \a) = r/l\0) is the lowest excited state 

For an oscillator with mass M, angular frequency fl, such that MVt/h = 1. 
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with one fermion [TTJ . The matrix element is evaluated by expanding af in diagonal 
fermionic operators. 

The magnetization profile m(l) can be expressed as a determinant [37] 



m(i) 



Hi G\ t i G12 ••• G\,i-i 

H2 G2,l G*2,2 ••• G2J-I 

Hi Gii G12 ■■■ Gii^i 



(3.14) 



where 



(3.15) 



n=0 



The rescaled magnetization profiles are shown in figure [2] where a good data collapse 
is obtained. The slopes are in good agreement with the values following from (12.141) and 
(12. 16ft : x ms — x m = 3/8 for z <C £, Luux m = 1/8 when to — 1 and 1/4 when u = 2 for 
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Figure 3. Finite-size scaling function of the surface magnetization for u> — 1 (left) and 
u) = 2 (right). The lines indicate the limiting behaviours given in (|3.18j) and (|3.19p . 



When 1 = 1, the determinant m(l) reduces to the normalized eigenvector <j)o(l) 
which can be evaluated for arbitrary values of the couplings [35l [36] . It gives the surface 
magnetization: 



0o(l) 



1 



L-l fc 

£11*? 

fe=i i=i 



-1/2 



In the unperturbed system this leads to m s ~ L 1 / 2 and 



(3.16) 

x m s = 1/2 as indicated in (13. 3p . 
In the scaling limit [L — > 00, £ — > 00, u(z) = Az/i) the surface magnetization is 
given by <f>o(u) at u(l) = so that, according to (13. lip and (13.121) . when to = 1: 

1)-. ,3,r) 
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Figure 4. Scaling function of the surface magnetization, reduced to its exponential 
part. In the semi-logarithmic scale, the decay is linear in |(;|L a ' +1 in agreement with 
the result of (I3.19P for g < (solid lines). 



The general case is discussed in the appendix where the finite-size behaviour of the 
surface magnetization in (13 . 16[) is studied in different limiting cases. When the scaling 
variable u(L) = \g\^'^ u+l ^ L <C 1 the leading behaviour is that of a finite-size critical 
system with 



m s 



L -i/2 



1 + sign(#)- v 



(oo + l)(w + 2) 

in agreement with (12.101) and (12.111) . 

When u(L) 3> 1 and L such that \g\L u = u u+i y L <C 1, one obtains 

\Of, , 4- 1^lV[2(w+l)] 



(3.18) 



m s ~ V 2 bl^ exp J +1 j , p < . (3.19) 

The result for m s in the ordered phase agrees with (13.171) when to — 1 and its in- 
dependence is the one expected from (12.121) . 

Numerical data for the finite-size scaling function of the surface magnetization are 
compared to analytical results in figure [31 The linear decay for L <C I has a slope 
equal to — x ms = —1/2 which follows from ( 12. lip and (13. 18ft . The constant values for 
L ^> £, given by (I3.19P for g > 0, are respectively ln[A m (l)] = 0.06039. . . for uj — 1 
and In [Am (2)] = —0.01098 . . . for u — 2. The exponential decay for g < and L ^> £ is 
compared to the asymptotic analytical result of equation ( I3.20p in figure [U 

3.4- Energy density profile and surface energy density 

The scaling behaviour of the energy density can be studied on the off-diagonal matrix 
element |(e|af|0)| where |e) = rj\ril\0) is the lowest excited state with two fermionic 
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-4-3-2-1 1 2 -4-3-2-101234 



ln(\g\ z) ln(\g\ Z) 

Figure 5. Scaling function of the energy density profile in the disordered phase. For 
uj = 1 (left) the solid line shows the profile obtained in (|3.21[) in the scaling limit. For 
uj = 2 (right) the solid line indicates the slope expected from (|2.14p when z <C I and 
the inset shows the exponential decay in the variable \g\z 3 . 




-4-3-2-10 1 2 -4-3-2-101234 

ln(g 1/2 z) ln(g" 3 z) 

Figure 6. Scaling function of the energy density profile in the ordered phase. For 
uj = 1 (left) the solid line shows the profile obtained in (|3.22[) in the scaling limit. For 
uj = 2 (right) the solid line indicates the slope expected from (|2.14[) when z <C ( and 
the inset shows the exponential decay in the variable \g\z 3 . The finite-size anomalies 
near the second surface at z — L are discussed in the text. 



excitations (see p2| for details). The expansion of of in terms of diagonal fermions 
leads to [381 



e{l) = \Ml)Ml)-Ml)<t>i(l)\- (3-20) 
When uj = 1, in the scaling limit with u = u(z) = \g\ l l 2 z, equations (I3.13p . (13.141) 
and (I3.15P can be used to calculate the energy profile which is given by 

eO) = \X3(u)xo(u) - Xi(u)x2(u)\ = A~\g\z I 1 + — =- e" 191 ^ , 
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when the system is disordered and by 

e(z) = \ X i(u)x (u)\ = Atgze-^ z \ At = 2^|, g > , (3.22) 

in the ordered phase. 

The numerical data for the scaled energy profiles are shown in figure [5] for g < and 
figure [6] for g > 0. A goog data collapse is obtained except for a finite-size anomaly near 
L when g > which is sent to infinity in the scaling limit. When to — 1 the convergence 
to the exact profiles obtained in the scaling limit is rapid. In the log-log scale the slope 
of the linear growth for z <C £ is equal to expected from (12.141) and 

obtained analytically in (13.211) and (13.221) . For z ^> £ the decay is exponential in the 
variable \g\z 2 . The insets show that for to = 2 the scaled profiles display a similar decay 
in the variable \g\z 3 . More generally one excepts an exponential decay in the variable 
\g\^ x/v according to (1231) . 
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Figure 7. Finite-size scaling function of the surface energy density for cu = 1 (left) 
and w = 2 (right). The lines indicate the slope -2 which follows from (|2.1ip for L <C ^ 
and the constant asymptotic values for L»l 



The surface energy density e s is given by equation (I3.20p at I — 1. Its finite-size 
behaviour is shown in figured When L <^ £ the slope is equal to expected 
from (12.111) . When L ^> £ the scaling function tends to a constant value which is given 
by g~ 1 e{z) at z — 1 in the limit g — > 0. For = 1 one obtains the constants v4~ in (13.211) 
and A~t in (I3.22p . In the logarithmic scale of figure [7] this leads to hiv4~ = —0.19112 . . . 
and lnA+ = 0.46736 . . .. 



4. Conclusion 



We have presented a theory for the scaling behaviour of physical profiles and finite-size 
behaviour in a semi-infinite [z > 0) critical system with free boundary condition in 
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the presence of a space-dependent deviation from the critical coupling. The coupling is 
critical at the surface z = and deviates into the bulk as a power law, A(z) = gz 1 ^ . The 
system is ordered (disordered) when the amplitude g > (g < 0). The perturbation 
introduces a new length scale, which is given by £ ~ anc l depends on the 

correlation length exponent v of the unperturbed system for a thermal-like perturbation. 
This length scale appears as a supplementary relevant size which must be considered in 
the scaling theory. 

The scaling arguments have been confirmed by an exact solution of the semi-infinite 
Ising quantum chain with free boundary condition in the case of a linear variation of the 
transverse field and checked numerically for the quadratic case. In the linear case, in 
the scaling limit L — > oo and g — > with Lg fixed, the excitation spectrum of the Ising 
chain is exactly mapped onto an harmonic oscillator problem with either Neumann or 
Dirichlet boundary conditions at the origin. Thus the magnetization and energy density 
profiles are given in terms of the harmonic oscillator eigenfunctions. 

In this study we have considered the effect of gradient perturbations at the ordinary 
surface transition. Other types of surface transitions, such as special and extraordinary 
transitions, can be considered as well. We are currently investigating these issues. 



Appendix 

The surface magnetization in (13.161) can be written as: 

L-l k 

S- l '\ S=l + J2Pk, lnft = X>(l-0n- (A - 1} 



L-l 

k=l 1=1 

When gL w <C 1 one obtains 



Jl_ nk^ 1 

(A.2) 

i=i u 

for k ^> 1. The sum in S can be replaced by an integral when L> 1, leading to: 

S ~ £ dk Pl ~ £ dk exp (-^y) ' (A ' 3) 
When gL w+1 1 the exponential can be expanded and integrated term-by-term 

s ^-^f d "" +1 -4-OTW' (A ' 4) 
m - i - 1/2 ( 1+ (^£T2))' < A ' 5) 

which gives equation (13.181) when the scaling variable u(L) is explicited. 

When gL w+1 ^> 1 and g > 0, neglecting an exponentially small correction, one may 
rewrite (1A.3D as 

dfcexpf-^-) ^JSL^I Ui (A . 6) 

Jo V W + 1 J [ 2 g(u + l)"] 1 '^ [2g{u + l)-] X /^ +1 ) 



and 
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which leads to the first expression of m s in (13.191) . 

Finally let us consider the case where gL w+1 ^> 1 and g < 0. Then in equation 
( 1A.3I) the argument of the exponential can be replaced by its first-order expansion near 
L, where the integrand takes its maximum value on the interval of integration 

f(k)~f{L) + {k-L)f(L), f(L) = ^^-, f'{L)=2\g\L u , (A.7) 

oj + 1 

and 

' L p/W 



S^e^ I dke (k - L ^' (L ^^—, (A.8) 



o 



where an exponentially small contribution from k = has been neglected. Equations 
( IA.3I) . (1A.7I) and (1A.8I) lead to the second expression of m s in ( 13. 19ft . 
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